00 



X 



General Non-minimal Kinetic coupling to gravity 

L.N. Granda*and W. Cardona^ 

Departamento de Fisica, Universidad del Valle 



o 

J5 Abstract 



We study a new model of scalar field with a general non-minimal kinetic cou- 



^^ pling to itself and to the curvature, as a source of dark energy, and analyze the 

' (— i ' cosmological dynamics of this model and the issue of accelerated expansion. A 

I wide variety of scalar fields and potentials giving rise to power-law expansion 
Q-i 

O have been found. The dynamical equation of state is studied for the two cases, 

'— ' without and with free kinetic term . In the first case, a behavior very close to 

CN that of the cosmological constant was found. In the second case, a solution was 
> 

\^ found, which match the current phenomenology of the dark energy. The model 

[^^ shows a rich variety of dynamical scenarios. 
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1 Introduction 

The astrophysical data from distant la supernovae observations [1], [2], [3], cosmic 
microwave background anisotropy [1], and large scale galaxy surveys [5], [6], [7], all 
indicate that the current Universe is not only expanding, it is accelerating due to 
some kind of negative-pressure form of matter known as dark energy ([8], [9], p^). 
The combined analysis of cosmological observations also suggests that the universe is 
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spatially fiat, and consists of about ~ 1/3 of dark matter, and ~ 2/3 of homogeneously 
distributed dark energy with negative pressure. This dark energy may consist of 
cosmological constant, conventionally associated with the energy of the vacuum or 
alternatively, could came from a dynamical varying scalar field at late times which 
also account for the missing energy density in the universe. That kind of scalar fields 
are allowed from several theories in particle physics and in multidimensional gravity, 
like Kaluza-Klein theory. String Theory[llj or Supergravity, in which the scalar field 
appears in a natural way. The cosmological solutions using scalar fields have been 
a subject of intensive study in the last years in the usual model called quintessence 
fn\ . [13], [E], [15]) [IS]) [H]) [IB]) where the stability of scaling solutions was studied 
(i.e. solutions where the scalar energy density scale in the same way like a barotropic 
fluid) for a field evolving in accordance with an exponential potential and a power 
law potential, which provide a late time inflation. More exotic approaches to the 
problem of dark energy using scalar fields are related with K-essence models, based on 
scalar field with non-standard kinetic term [19], [20]; string theory fundamental scalars 
known as tachyons pT|; scalar field with negative kinetic energy, which provides 
a solution known as phantom dark energy [22] (see [8] for a review). The non- 
minimal coupling between the quintessence field and curvature, as an explanation of 
the accelerated expansion, has been considered among others, in refs. [23] [21], [25]. 
This coupling has also been studied in the context of the inflationary cosmology [26] , 
|27j . [2H]) [22]) [3D]- An inflationary model including the DBI term and non-minimally 
coupled scalar has been proposed in [311 [32] • 

Considering the possibility of having a situation where the accelerated expansion 
arises out of modifications to the kinetic energy of the scalar field, in this paper 
we consider an explicit coupling between the scalar field, the kinetic term and the 
curvature, as a source of dark energy, and analyze the role of this new coupling in 
an evolution scenario with late-time accelerated expansion. The basic motivation for 
studying such theories is related with the fact that they appear as low energy limit 
of several higher dimensional theories, e.g. superstring theory [H], and provide a 
possible approach to quantum gravity from a perturbative point of view. They appear 
as part of the Weyl anomaly in A^ = 4 conformal supergravity [33] IM]. A model with 



non-miniinal derivative couplings was proposed in [35], [36], [37] in the context of 
inflationary cosmology, and recently, non-minimal derivative coupling of the Higgs 
field was considered in [38], also as inflationary model. In [39j a derivative coupling 
to Ricci tensor has been considered to study cosmological restrictions on the coupling 
parameter, and the role of this coupling during inflation. Some asymptotical solutions 
for a non-minimal kinetic coupling to scalar and Ricci curvatures were found in [ID] , 
and quintessence and phantom cosmological scenarios with non-minimal derivative 
coupling have been studied in [H]. Non-minimal coupling of arbitrary function f{R) 
with matter Lagrangian (including the kinetic scalar term) has been introduced in 
[12], [13]. Such a model was proposed to describe the dark energy, late-time universe 
acceleration. When f{R) represents the power law function, such a model (which 
maybe considered as string-inspired theory) maybe proposed for dynamical resolution 
of cosmological constant problem as it has been suggested in 01], [15]. In the present 
work we consider this function f{R) as linear in R but we generalize the model 
permitting extra i?^,^ coupling with kinetic-like scalar term. In addition, we keep the 
free kinetic term for the scalar field. 

We propose a model of dark energy in which the scalar-field kinetic term is non- 
minimally coupled to itself and to the scalar and Ricci curvatures. We consider a 
general self coupling F{(j)), which allows to find a wide variety of potentials giving 
rise to accelerated expansion. Some solutions for the product F(0)0^ are proposed, 
and the scalar field and potential are found. We have considered power-law solutions, 
and solutions giving rise to dynamically varying equation of state. In the last section 
we present some conclusions. 



2 Field Equations 

Let us start with the following action, which generalizes the model proposed in 
S = J d'x^^g [^R - \d^^d^^ - \iR (F(0)9^09» - 



\r^, {F{<p)d^<pd''<P) - V{<P) 



(2.1) 



where Sm is the dark matter action which describes a fluid with barotropic equation 
of state. The dimensionahty of the couphng constants S, and rj depends on the type 



of function F{(f)). Taking the variation of action 2.1 with respect to the metric, we 
obtain a general expression of the form 

1 



^/lu ^Ofiu^ — ^ [-'■ fiu ~^ -'m'^J 



(2.2) 



where k^ = SttG, T™ is the usual energy-momentum tensor for matter component, 
and the tensor T^j, represents the variation of the terms which depend on the scalar 
field d) and can be written as 



T, 
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rrifj) _|_ T^C _|_ rpV 

ii.u ~i~ itu ~r -^ 1 1.1 J 



where T^^, T^ui '^uu correspond to the variations of the minimally coupled terms, the 
^ and the r] couplings respectively. Due to the interacting terms this expressions are 
defined in the Jordan frame and do not correspond to the energy-momentum tensors 
as defined in the Einstein frame. Those variations are given by 



7;t = V^0V,0 - -g^,Vx<pV'<P - g^,V{<p) 



- IjLV 
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Tl =e R,u - ^9,uR (F(0) Va0V^</)) + g^uVxV^ {F{<j>)V,<j>V^ 



-iy^Nu + V,V^) (F(0)Va0V"0) + R (F(0) V^0V, 



(2.4) 



(2.5) 



T;, =r/ F(0) (i?^AVVV.0 + i?.AVVV, 



-g^uRx, (F(0)VVV 
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- 2 i^>^^^ (F(0)V>V.0) + VaV. (F(0)V>V^0)) 

+ ^VaV" (F(0)V^0V,0) + ^(7^.VaV-, (F(0)VVV^0; 
Variating with respect to the scalar field gives rise to the equation of motion 

^ -.d, [v^ (ei?F(0)9^0 + r^R^^''F{<P)d,cP + 9^0)] + ^+ 



(2.6) 
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dF 



(2.7) 



md, 



+ 7]Ruudyd' 







Assuming the spatially-flat FriedmannRobertson Walker (FRW) metric, 

ds^ = -de + a{tf [dr^ + r^n'') 



{21 



From Eqs. (2.4 2.6) and using the metric (2.8) we can write the (00) and (11) com- 



ponents of the Eq. (2.2) (with the Hubble parameter if, and for homogeneous time- 



depending scalar field) as follows 



H' 



K 



2^1 V 



+ V{4>) + 9^H^F{(t)),p^ + 3(2^ + r])HF 



3 dF ■ 

- 3(2^ + v)HF{<p)<p<p - 2 (2e + V)H^<P' 



(2.9) 



and 



-2H - 3H^ =K^ 



rl 



02 - V(0) + 3(e + v)H'Fi<P)^' + 2(e + v)HFicj>)^' 



+ 4(e + v)HF 

+ (2^ + r^)(F{<P)4>' + F 



+ 2{^ + 7])H^<j)^ 



(2.10) 
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5dF 



+ 



Id^F :, 



2 rf0 " " 2d(f)^' 
where we have assumed scalar field dominance (i.e. TJ^ = 0). The equation of motion 



for the scalar field (2.7) takes the form 



dV ... . . 

+ 3H(t) +-T- + 3(2^ + 7])HF{(t))(t) + 3(14^ + 5r])HHF{4>)4> 



+ -i2^ + v)H(2F 



dF •, 



-(4^ + r])H^2F{4>)4> + 



dF ■.. 



(2.11) 



+ 9(4^ + r])H^F{(f))^ = 

where the first three terms correspond to the minimally coupled field. In what follows 
we study the cosmological consecuences of this equations, under some conditions that 
simplify the search for solutions, but nevertheless show the richness of the cosmological 
dynamics of the present model. 



The Eqs. (2.9 2.11) significantly simplify under the restriction on ^ and r] given by 

r/ + 2^ = (2.12) 



In this case the field equations (2.9 2.11 ) contain only second derivatives of the metric 
and the scalar field, avoiding problems with higher order derivatives [361 HO]- The 



modified Friedmann equations (2.9) and (2.10) take the form 

,2 



H' 



K 



1 • 



0' + V(0) + 9^H'F{(f))(f)' 



(2.13) 



and 
-2H-3H^ = K^ 



1 •, 



dF -. 



- V{<P) - ^ {3H' + 2H) F{<f))^' - 2iH f 2F(0)00 + ^ -^^ 

(2.14) 
where we have replaced rj = —2^, and the equation of motion reduces to 



+ 3/f0+ ^ + 3^^' (2F(0)0 + ^0M + 18^H^F{(P)^+ 12^HHF{<P)^ = (2.15) 

Next we try to find cosmological solutions, giving rise to to accelerated expansion 
and acceptable behavior of the equation of state parameter (EoS). In the point A of 



the next section we consider the general equations (2.9) and (2.11) in the special case 



when the kinetic coupling dominates over the free kinetic term. But in the rest of the 



paper we will study the cosmological implications of the Eqs. (2.13) and (2.15) 



First note that, independently of F{(f)), if we consider the asymptotic solution = 
00 = const., then from Eqs. (2.9[) and (2.11) it follows that V = Vq = const and 



H = Hq = Ka/Vo/3 and we get an asymptotic de Sitter behavior. 



3 Power-law solutions without free kinetic term 

It is of interest to derive a scalar coupling function F{(f)) that gives rise to a power-law 
expansion (a oc t^), as it's well known that this kind of expansion is characteristic of 
the evolution at the early stages of radiation and matter dominance, and also may 
respond for accelerated expansion. 



Let's consider the particular case of the model (2.1 ) without free kinetic term. In this 



case we are in the frames of above models (|l2l |13], HU HH]) which correspond to pure 
dark energy scalar sector. This also applies if the the slow-roll condition cfp' « 1^(0) 
is considered, which can take place in the contexts of dark energy or infiationary 
cosmology. Let us consider the effects of this new kinetic coupling in the cosmological 
dynamics, in the case of scalar field dominance (we further neglect any background 
radiation or matter contribution). 



A. Derivative coupling with t] + 2^ ^ 



Let us then, begin with the simple case of strictly non-minimal kinetic coupling, 
without potential term. Introducing the new variable x = '^^-^(0)0^5 and considering 



first the general Eqs. (2.9 2.11) with arbitrary C, and t] and without free kinetic and 



potential terms, the Friedmann equation and the equation of motion from Eqs. (2.9) 



and (2.11) take the form 



1 



H'' = 3^H\ + (2^ + r])Hx " §(2^ + v)Hx 



{A^ + r])H^ + {2^ + r])H 



X + 2 



3(4e + v)H^ + (14^ + 5v)HH + (2^ + v)H 



(3.1) 

x = o 

(3.2) 



where we have multiplied the Eq. (2.11 ) by 0. Let's propose now a power law solution 



H = p/t, and replace in Eq. (3.2) to obtain 



[(4^ + r])p - (2^ + v)]x + '2 [3(4^ + v)p' - (14^ + 5v)P + 2(2^ + r])]x = (3.3) 



The solution to this equation is of the form 

3(4^ + v)P^ - (14^ + 5r])p + 2(2^ + v) 



x = t' 



a 



(3.4) 



2^ + r]-{A^ + 7])p 
where we considered the condition xo = 1 at to = 1- Replacing this solution in Eq. 



(3.1), we obtain the following restrictions on the parameters 

3(4^ + r])P^ - (14e + 5r])p + 2(2^ + r^) = 0, 
P 



2^ + V 



(3.5) 



3^-1 
giving rise to the simple restriction for accelerated expansion {p > 1): ^ > 1/3, ^ — r/ < 

1 or ^ < 1/3,^ — rj > 1. Replacing p from the first condition (which is equivalent to 

a = 0) in the second one, gives an equation for ^ and r], which has the three following 

solutions 

77 = -2/3, r7 = -2e, r/ = -e-l (3.6) 

the solution rj = —2/3 can not give rise to accelerated expansion, the second solution 
gives p = 0, and the third solution gives accelerated expansion on the line rj = —^ — 1 



in the interval < ^ < 1/3. The scalar field can be obtained from the Eq. (3.4) for 
a = 0, as follows 

K I ^F{(l))d(t) = ±t (3.7) 

for a given scalar function F{(j)). Thus for -F(0) = 1/0^, as considered in |16], and 
taking a = 0, the scalar field is obtained as (f) =oc e^*'''^. The above simplifications 
allowed us to check that the power-law solution takes place even without the restric- 



tion (2.12) 



B. Derivative coupling with 77 + 2,^ = 



From now on, we will assume that the relation (2.12) holds. Turning to the equations 



(2.13)-(2.15), the Eq. (2.13) reduces to 



3^K^F{4>)^' 



(3.8) 



which can be solved in terms of the integral 



y^FW)d(l> 



3^K 



(3.9) 



and the Eq. (2.15) (without free kinetic term and potential) takes the simple form 



3H' + 2H = 



(3.10) 



and therefore H = 2/(3t), which gives a power-law corresponding to pressureless 



matter dominance. Note that substituting back to Lagrangian the solution (3.8), 



(3.9), we get kind of non-covariant modified gravity (for general review , see ( |47j ) ) . 



From Eq. (2.13), the effective gravitational coupling can be expressed as 



nljf = 8nG,ff = k\1 - 3^K'^'F{<j))r' 



(3.11) 



by other hand, taking the derivative of the above expression, and using the equations 
for X = k20^-F(0) given in (3.4) {a = 0) and (3.8), it follows that Geff/Ge/f = 0, 



satisfying the restrictions on the time variation of the gravitational coupling |48 
C. Derivative coupling with potential 



Let's now consider the scalar field potential. We will look for the shape of the 
potential V{(j)) corresponding to the asymptotic behavior giving rise to accelerated 



expansion. In this case, adopting the same variable x = i^'^F{'P)'P'^ ■, the Eq. (2.13) 
takes the form 



H 



^_K^ v{4>) 



3 1 - sex 



and multiplying by k (p the Eq. of motion (2.15), it's obtained 



AV 



K 



2^^ , ^^H'^^^i^^H^^^i2^HHx = Q 



dt 



dt 



(3.12) 



(3.13) 



where the first two terms of Eq. (2.15) have been dropped due to the absence of the 



free kinetic term. Taking the derivative of the potential in Eq. (3.12), and replacing 



in Eq. (3.13), we obtain 



ei/'^ + iHHx - S^H^x -HH = 



looking for power-law solutions, we replace H = p/t in Eq. (3.14), obtaining 



et^-^(l + 3p)x+l = 



(3.14) 



(3.15) 



considering the particular solution x 



c(i]^ ) , gives the equation for the scalar field 



K^F 



(3.16) 



e(l + 3p) 
which gives the scalar field for a given function F{(f)). For instance, if we take F{(f)) 

1/0^ as in |in], then 

t 



5oexp 



/^Ve(l + 3p) 



(3.17) 



And the time dependence of the scalar field potential is obtained from Eq. (3.12) by 

replacing H and x 

3p2 / 3j5 - 2 



note that despite the fact that the time dependence of the potential is known, the 



shape of the potential depends on the coupling -F(0) through Eq. (3.16). In fact, for 



-W-^(i|:=^)i^ (3.«) 



F(0) = 1/02 as in 



Note that we have not restriction on p, so for p > 1 it follows the accelerated expan- 
sion. 
Another interesting coupling is given by F{(f)) = e'^^'^/(f)Q with constant A. Then from 



(3.16), the scalar field is obtained as follows 



H^vmm^ 



and the potential is given by 



_ 3AVy(3p-2) 

this field and potential generates power-law expansion in the quintessence model 
(|49tlS|). We can also obtain a potential of the inverse power-law type, giving rise to 
power-law expansion. If we consider the coupling F{(f)) = M~2~"(/)", then 

f n + 2 \^ 2 , ^ 

\2k.JW+W)) 

and the potential is of the form 

_ / 3p2(3p-2)(n + 2)2 \ M"+2 
^" V 4e«:4(l + 3p)2 J 0"+2 ^^-^^^ 

Note that the case n = reproduces the tachyon potential giving rise to power-law 

expansion (see [I^ [5]) 



The general solution of Eq. (3.15) is given by 

X = n'Fm' = e^'^ + ^^^^ (3.24) 

which after integration with respect to time, gives the expression for the scalar field 
as follows 

^ 3(i + p)(m + 3p))^/ 



1/2 



tm+m 



1/2 



H r — Hypergeometric2Fl 

3^(1 +p) 



11. 1 



l + ^^^,-e(l + 3p)t 



l+3p 



l + 3p 2' l + 3p 

(3.25) 
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The time dependence of the potential, as follows from (3.24) and (3.12) (with H = p/t) 
is given by 



V{t) 



3p^ f3p-2 

3p+l' 



r 



3^^3p-l 



(3.26) 



which gives the shape of the potential V{(j)) through Eq.(3.25). For such configuration 
of the scalar field and potential, and despite the fact that the potential can not be 
solved explicitly in terms of the field, it is possible to have power-law expansion with 
p > 1. The asymptotic behavior of this potential at t — )■ and t — )■ oo is the same 
as discussed in [16] for F{(j)) = 1/0^. The only difference is that the shape of the 



potential depends on -F(0) through Eq. (3.25). 

One can use the freedom in choosing the coupling function -F(0) in order to simplify 

the integration in Eq. (3.24). Making n = 1 and proposing a scalar field of the form 



t, the Eq. (3.24) becomes 



l1+3p 



+ 



e(l + 3p) 



1/2 



dt 



which gives a solution 



F(0) 



a+3p 



+ 



1 
e(l + 3p) 



(3.27) 



(3.28) 



And the corresponding potential, as seen from (3.12) and (3.24), can be expressed 
explicitly in terms of the scalar field as 

2/'3p-2 1 



y(0) = 3p' 



3^0 



3p-l 



(3.29) 



So this coupling function gives rise to power-law expansion a (x t^, without any 
restriction on positive values of p. 



For the particular solution Eq. (3.16), the restriction on the time variation of the 



Newtonian coupling is satisfied, as discussed at the end of the point B. For the 
general solution Eq. (3.24) we can follow the same discussion given in [IH] (see Eq. 
(3.24) in US]). 
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4 Dynamically varying equation of state without 
free kinetic term 

In this case we propose an asymptotical solution for the scalar field at late times, and 
try to reconstruct the Hubble parameter and the potential, giving rise to dynamically 
changing EOS parameter, evolving according to the current observations. Let us write 



the equations (2.13) and (2.15) in terms of the variable x = log a. The Eq. (2.13) 



without free kinetic term becomes 



H'='^ (r(0) + 9^H'<p"Fi<P)) 



and the Eq. of motion (2.15), after multiplying by takes the form 
dV 



dx 
where ch' 



+ 3^HU2<p'<f)"F{<p) + <p''^ + 



dx 



"F{<p))+m'^<P"F{<p)=0 



[A.l] 



(4.2) 



dcp/dx, and we used J^ 



dx 



Resolving the Eq. (4.1) with respect to H^ one obtains 

1 



H' 



(l ± V^l - A^K,^V(j)'^F{(j)) 



(4.3) 



In this case, in order to keep H^ positive, we should take care of the sing of the 
combination ^F (assuming that V is positive). So, if ^F > 0, we can take both signs 
of the root, and for C,F < we take the (— ) root. In any case, the choice of the 



sign does'n infiuence the Eq. (|4.2|), as V in Eq. (|4.3|) is under the square root. The 

(4.4) 



simplest case of Eqs. (4.2) and (4.3) is obtained by considering (see 

0'^F(0) = const. = c? 



here we assume that -F(0) > 0. Replacing in Eq. (4.3) we get 

1 



E' 



«:^~. 



V[x), V{x) 



lia^K 



2^4 



1± Vl-4^«^/€^1/ 



(4.5) 



setting K^ = 1 and replacing Eq. (4.5) in Eq. (4.2), we obtain the following equation 
for the new potential V[x) 



ia^V{x)-lY^-2ia^V\x) = ^ 



(4.6) 
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Integrating this equation we obtain the solution for V as follows 



V{x) 



^a^ProductLog 



-2x-l/i^a2Vo) 



(4.7) 



using this solution in Eq. (4.5) (with k^ = 1), and changing to the redshift variable 



z, the Hubble parameter is obtained as follows 



H\z) 



S^a^ProductLog 



1 



-y(^-'y")(i + zf 



{A.l 



Using the properties of the ProductLog function, and the sign of the coupling constant 
^ (here Vq > 0, and without loss of generality we assumed F{(j)) > 0, as the kinetic 
coupling is given by the product ^F), we can follow the same arguments and results 
given in [36], which I resume here: For positive values of the constant ^ we can enter 
in a forbidden region for the arguments of the ProductLog function. In order to avoid 
those regions where the ProductLog function takes complex values, we consider ^ < 
which will give always H^ positive (see discussion in 



From the Eq. (4.5) and the expression for H^, it follows that the effective dark energy 
density given by p^ff = SH'^/k,^ = V tends to zero at high redshifts, validating the 
assumption that the dark energy was negligible at early times, when the matter 
contribution was relevant, but presents a Big Rip singularity at the future, revealing 
phantom behavior. The effective equation of state We/j = — 1 — ^-^^^r, from Eqs. 



(4.5) and (4.7) in terms of the redshift, can be written as 

2 



w. 



eff 



-1 



3 ( ProductLog 



(a^Vo 



e-V(5"^^o)(i + ;2)2l +ij 



(4.9) 



which varies very slowly with z for a given combination of ^,a,VQ. Fig.l shows 
the behavior of the effective equation of state for two different combinations of the 
parameters, represented as 7 = ^a^Vo. 
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Fig. 1 The dark energy effective equation of state parameter Weff versus redshift, for 
two combinations of the parameters 7 = —0.1 (top), and 7 = —0.5 (bottom), where 

7 = ^a'Vo 

For the product of parameters 7 = ^a^Vo = —0.1, the effective EoS parameter varies 
from Weff ~ —1.06 at the present {z = 0) to —1 at 2; — )■ 00, showing very slowly 
variation with time, and maintaining very close to —1. This behavior is very close to 
that of the cosmological constant. 

From Eq. ( |4.4 ), for a given coupling function we can find the scalar field. Thus for 
F(0) = 1/0^, solving (4.4), the scalar field depends on a; as = 0oe~°^. Assuming 
F(0) = (J)"', the scalar field becomes (p = -^^t2/("-+2) 



The time variation of the effective gravitational coupling can be found, if we use the 
fact that F((h)(h''^ = o? , which allows to write 



G 



eff 



G 



aHn 



a 



Q^a^K^H^ 



(4.10) 



eff to=i 1 — S^a^/t^ifg 

considering Hq < iJp, and due to the extremely small value of n^H^ = MpfH"^, is al- 
ways possible to find a suitable combination of ^a^ that satisfy the current constraints 

or^Glff/GeffM- 



5 Cosmological solutions with free kinetic term 

In this section we include the free kinetic term and analyze some important results, 
specially concerning the dynamical equation of state. 
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A. De Sitter Solution 



Multiplying the Eq. (2.15) by 0, and replacing the product F{(j))(j)'^ from Eq. 



(2.13), the Eq. (2.15) reduces to first order equation with respect to the variables 
if) = (f?, H and V , and can be written as 

H^ + (qH'^ -H\i) + 2H-^ - 2 (3H^ + h) V + 6^ (^H^ + 2H\ = (5.1) 

Using the fact that the variables ip and V are separated, we can limit the model to 
the class of potentials that satisfy the restriction 



H—- - [?,H^ + H]V + 3^r ( 3if ^ + 2H 







(5.2) 



dt \ / K^ 

then, the equation for the field ip significantly simplifies, but still giving an interesting 
dynamics as can be seen bellow 



H^+ {6H^-H)iP = 



(5.3) 



Assuming the de Sitter solution a(t) = aoe^°^, the solution to Eqs. (5.3) and (5.2) 
for the scalar field and potential, are given by 



ioexp(-3/Jot), V = Vo^ + ^ 

(h K^ 



(5.4) 



and replacing in Eq. (2.13), the following expression for the coupling function is 
obtained 



n<t>) 



1 



21^Hl 



1 



n^(j)^ 3p^ 



log^ 



+ 



ioM) ^ 1 



2 ■ 9iJoV^ 3fi;202 



(5.5) 



In this manner, the solutions given by (5.4) together with the coupling (5.5) give an 



exact de Sitter (inflationary) solution, which in the case of the standard nonminimaly 
coupled scalar field is obtained under the known slow-roll conditions. 
B. Dynamically varying equation of state 



Considering the full model with the restriction (2.12), we will use the freedom in 



choosing the coupling function in order to obtain a cosmological dynamics more close 
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to the one expected from observations. In terms of the variable x = log a and defining 
0'^, the Eq. (2.15) can be written as (after multiplying by 4>) 



the function 9{x 
1 d 



2dx 



[H^e) + ?,H^e + ^ + 9^H^^Fe + 3^H^^{Fe) + IS^H^FO = (5.6) 



dx 



dx 



dx 



From Eq. (2.13), changing to the variable x, we can write the product F0'^ = F6 as 
following 



F9 



1 



9 



V 



3^k'^H'^ 18^H^ %H^ 



(5.7) 



taking the derivative of Eq. (5.7) and replacing F6 and d{F6)/dx into Eq. (5.6), we 



arrive at the following equation involving 6, H and V 

,de ,/ , dH'^\ ^dV ( , dH'^\ 

dx \ dx J dx \ dx 



+ 12- 



H^ 



3H' + 



dH^ 

dx 



{5.t 







In this manner, we obtain a first order differential equation for the functions 6, H 
and F. In order to integrate this equation, and thanks to the fact that the functions 
9 and V are separated, we can impose a restriction on the scalar field potential given 
by the equation 



2^2^ -{qH^ + —\ V + 6^H' (3H' + —1=0 



dx 



dx 



K^ 



dx 



(5.9) 



which simplifies the Eq. (5.8): 



.-vo de 



2H' 



dx 



+ 12H^ + 



dJP 

dx 



= 



(5.10) 



where we have defined the scaled Hubble parameter H^ = H'^/Hq (note that this 
equivalent to divide the Eq. (5.8 by Hq). Let us propose an interesting solution for 



the Hubble parameter, which is important in describing the dark energy, and can fit 
the observations (see [SU]) 



H^ = fi^e""^ + nBe-^"", ^A > 0, fis > 



(5.11) 
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where the diinensionless parameters Qa and Qb satisfy the natural condition at the 
present time {x = 0) 

nA + nB = i (5.12) 

with this proposal for H^ we can find the functions 6{x) and V{x), that satisfy the 



Eqs. (5.9) and (5.10). The constants Qa, cx and /3 can be fixed by comparing the 



behavior of the effective equation of state with that expected from observational data. 
Replacing H^ in (5.10), after integration we obtain the following expression for Q 



e{x) = ^' = 00 



62 



|{a+/3-12)x 



(fi^e/3- + nee^ 



,1/2 



(5.13) 



where ^o is the integration constant. Replacing now H^ in (5.9), and rescaling the 
potential in the form V = k^V/Hq, it is found the following expression for the scaled 
(dimensionless) scalar field potential 



V{x) = Cie 



-Ua+l3-&)x 



1/2 



(n^e^" + ^Be"") ' + 6e 



(l + g^e^(- 



1/2 



VLA{a — 3) (2/3 — a + 6)Hypergeometric2Fl 
Vtsia + 6)(/3 - 3)e(°-^)'^Hypergeometric2Fl 



(a-l3)x 



(a + 6)(a-2/3-6) 
6 + a 1 a-2(/3 + 3) VLb 

a-2(/3 + 3) 1 4/3 - 3a + 6 _^ („_ 
2(a-/3) '2' 2(/3-a) ' fi^^ 



_2(/3-a)'2' 2(a-/3) 



+ 



(5.14) 
where the integration constant Ci is now dimensionless. Integrating the square root 



of Eq. (5.13), we find the scalar field as 



1 + 



QA\ 



^A J (fi^e/^^ + fise"^)^/^ 
Hypergeometric2F 1 



a -12 1 a -12 n 

4(a-/3)'4' ^4(a-/3)'~a^'' 



-B (a-/3)a; 



(5.15) 



Finally, the coupling function F{x), is obtained by replacing the Eqs. (5.11 5.14) in 



(5.7). 



Hence, starting from the model with the action (2.1) (with Sm = 0), as we show 
bellow, one can find an exact solution that not only gives the accelerated expansion, 
but also gives the correct transition from decelerating phase to accelerating phase (for 
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transition deceleration-acceleration in /(G)-gravity see |51j). 
Expressing H^ in terms of the redshift z, it takes the form 

writing the effective equation of state parameter (EoS) in the form 

^ 1 (1 + -2) d{H^) 



(5.16) 



w. 



eff 



(5.17) 



3 m dz 

we can adjust the constants a, /3 and VLa in order to obtain the appropriate behavior, 
compatible with the observational data. In Fig. 2 we plot the effective EoS for three 
different sets of the constants, and in the case corresponding to the set of parameters 
a = 3, /3 = —1.49 and VLa = 0.27 the (red) curve crosses the —1 divide at z ^ 0.1 
showing quintom behavior (for a recent review on observational and theoretical status 
of quintom models, see [52]). Note that all curves correctly describe the matter 
dominance at early times, the transition to accelerated expansion at Zt ~ 0.5 — 0.7 
(see [S3]), and the current accelerated expansion with the equation of state very close 
to -1. 




Fig. 2 The dark energy effective equation of state Weff versus redshift, for three sets 
of parameters: green-(a = 3, (3 = —1, Qa = 0.3), hlue-(a = 3, /3 = —1.49, VLa = 0.35J 

and red-(a = 3, /3 = -1.49, VLa = 0.27^ 
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In terms of x, the time variation of the gravitational couphng can be written as 



(-Jeff oqK d 



Geff 1 - S^k'^FHW dx 



{FH^d)H (5.18) 



replacing the numerical values of the parameters a, /3 and fi^, in (5.11), (5.13) and 



(5.14), we obtain the expressions for iJ, 9 and V depending on x and the constants 
Ci and 6*2 (where we replaced 6o = C2, which is a constant of dimension (lenght)'"^). 



Replacing this expressions in Eq. (5.7) and multiplying by H^, we get the product 
^K,'^F{x, Ci,C2)H{xY6{x, C2) , which is relevant for our analysis. After some algebra 
and numerical evaluation of the Hypergeometric2Fl function, we arrive at the result 



G. 



eff 



^*'^"^^) -H, (5.19) 



i„ l-3/(Ci,C2)' 

here we can use the constants Ci and C2 to guarantee that g{Ci, C2) ~ 10^^ or less, 
and f{Ci,C2) << 1 in order to accomplish the restrictions imposed by the current 
value and the time variation of the gravitational coupling. For instance, if we take the 



2 



set of parameters (a = 3, /3 = -l,fiA = 0.3), then g ^ -10.71 - 0.34Ci + 0.33k^C. 
and / ~ —3.42 — 0.1 ICi — 0.055fi;^C2 (note that C2 appears in the dimensionless 
combination k,^C2). Instead of inequalities, we can also satisfy the restrictions by 
solving the linear system g{Ci,C2) ~ and f{Ci,C2) ~ 0, which in this case gives 
Ci ^ -30.9 and C2 ^ 0.2. 



The solution 5.16 is very attractive, as is able to describe the current observational 
value of the EoS parameter, and the expected dynamical evolution of the EoS, show- 
ing not only the transition deceleration-acceleration, but also, for some values of the 
parameters the model shows quintom behavior. 

In Fig. 3 we have plotted the scalar field potentials corresponding to the set of param- 
eters given in Fig.2, making the same correspondence with the colors, and with the 
constants Ci and C2 satisfying the conditions g{Ci, C2) ~ and /(Ci, C2) ~ 0. Note 
that all the potentials present a very nice behavior, being negligible at high redshift, 
where we expect dark matter dominance, and becoming important as the redshift 
decreases, where we expect the main contribution of dark energy. 
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Fig. 3 The dark energy scalar field potentials in terms of x = log a (a^ = 1), for 
the three sets of parameters: green-(a = 3,(3 = —1,Qa = 0.3, Ci = — 30.9J, blue- 
fa = 3, (3 = -1.49, ^A = 0.35, Ci = -1479; and red-(a = 3,/3 = -1.49, fi^ = 
0.27, Ci = -1892; 

6 Discussion 

We have proposed a new model of non-miniinally kinetic coupled scalar field, where 
the kinetic term is not only coupled to itself through the function -F(0), but to 
the curvature, giving rise to interesting cosmological consecuences. An interesting 
characteristic of the model is that the power-law behavior depends on the dynamics 
generated by the the nonminimal kinetic term as a whole, and the coupling function 
F{(f)) is important in defining the time dependence of the scalar field and the shape 
of the potential. Even in the absence of a potential, the system has power-law so- 
lutions. First we considered the strictly non-minimal kinetic coupling, without the 
restriction on the coupling constants ^ and rj (i.e. rj ^ —2,^) and have found late time 
power-law solutions, leading to accelerated expansion (see Eq. (3.4 3.6[ )). With the 
restriction (r^ -|- 2,^ = 0), the strictly non-minimal kinetic coupling produces the effect 



of pressureless dark matter (3.10), which is the behavior of the matter at early stage 



of the evolution. In the presence of a potential V^(0), the model presents power-law 
solutions giving rise to accelerated expansion, as showed in the point C of section 3, 
for different couplings and potentials. In particular, the potentials for power-law so- 
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lutions of miniinally coupled quintessence and tachyon models are obtained in (3.21) 



and (3.23) (see [IHIIH])- Additionally, the potential (3.26) may be interpreted in two 



asymptotical cases, corresponding to early and late time behavior. 

Looking for solutions in the presence of potential, but giving rise to dynamical effec- 
tive EoS (not of the power-law type), we proposed a restriction for the scalar field and 
the coupling, of the form 0'^F(0) = const. = a"^, and found the potential presented 



in Eq. (4.7). According to this result, the model presents phantom behavior, but 
with an EoS very close to that of the cosmological constant, as seen from Fig. 1. 
Note from Eq. (4.5) and Fig.l, that when the combination 7 = —^a'^Vo — )■ —0, the 



effective EoS parameter We/f — ?■ — 1, behaving very close to the cosmological constant. 
The effective EoS parameter varies from —1 in the far past ^ — )■ 00 to —5/3 in the 
future [z = —1), but is worth to note that Eq. (4.5) allows any current value of We//o 



in the limits —5/3 < Weffo < "~1; and this value is closer to —1 for small values of 
7. For instance, for 7 = —0.1, the effective EoS parameter changes from —1.06 to 
— 1 in a practically infinite time, giving very exact description of the cosmological 
constant. In forthcoming work we are considering the matter contribution in order 
to fix the parameters of the model, according to the observational data. Note that 
in models with non-minimal coupling to f(R), we can consider also f{G) (where G 
is Gauss-Bonnet combination) as function multiplied to kinetic scalar term. Then it 



would of interest to study generalization of model (2.1) where instead of R we use 
function G (see ([S3])). 



By imposing the restriction given by Eq. (5.9) on the scalar field potential, we 
have found a class of solutions which describe an appropriate cosmological evolution 
as shown in Fig. 2, showing the transition from the deceleration to acceleration phase 
at zt ~ 0.5 — 0.7, and also crossing the phantom divide at z ~ 0.1 (the red curve). 
All the presented solutions respect the restrictions imposed by the time variation of 
the Newtonian coupling. 

In conclusion, the scalar field model with non-minimal kinetic coupling to itself and 
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to the curvature, provides a new framework to describe the actual phenomenology 
of the dark energy. The derivative couplings provide an effective scalar density and 
pressure which might play an important role in the explanation of the dark energy or 
cosmological constant. This makes the present model very attractive for future study 
of the dark energy problem. 
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